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1 Introduction 

Heavy quark effective theory (HQET) [|IJ provides a framework to study the heavy hadron 
spectra and transition amplitude with the systematic expansion in terms of 1/mq. Using 
the light cone photon wave function (PWF), radiative decay processes like B — > Ivy, 
B — ► p7 have been studied p], [| |J ^ || with QCD sum rules (QSR) @. Recently similar 
approach was employed to analyze the couplings of pions with heavy hadrons |§| . In this 
work we will study the radiative decays of heavy baryons and lowest three doublets of 
heavy mesons using the light cone QSR (LCQSR) || in the leading order of HQET. With 
LCQSR the continuum and excited states contribution is subtracted more cleanly, which is 
in contrast with the analysis of meson radiative decays using the external field method in 
QSR |H|. The LCQSRs for the radiative decays of heavy baryons are presented in section 
||. Section [| discusses radiative decays of lowest three doublets of heavy mesons. The 
following section is a discussion of the parameters and the photon wave functions. The last 
section is the numerical analysis and a short summary. 

2 Radiative decays of heavy baryons 

We first introduce the interpolating currents for the heavy baryons: 

Va(x) = e abc [u aT (x)C l5 d b (x)}h c v (x) , (1) 
rfr{x) = e^lu^^C^d'ixM^Kix) , (2) 



i£h..(s) = e abc [u aT (x)C lu u b (x)}(- 9 r + \lW)K{x) , (3) 

where a, b, c is the color index, u(x), d(x), h v (x) is the up, down and heavy quark fields, T 
denotes the transpose, C is the charge conjugate matrix, g^ v = g^ v — v^v u , 7^ = 7 M — vv M , 
and v M is the velocity of the heavy hadron. 

The overlap amplititudes of the interpolating currents with the heavy baryons is defined 

as: 

(0Na|A)=/a«a, (4) 
(0|7teiE) = / E u Ej (5) 

(0|^|£*> = ^|<», (6) 

where w s * is the Rarita-Schwinger spinor in HQET. In the leading order of HQET, / E = / E * 

The coupling constants 77^ are defined through the following amplitudes: 

M(S C -> A c 7) = ier] 1 u Ac (T' J ' 1 'q fl e 1/ ux c , (7) 
M(£* -> A c7 ) = ^er/ 2 e^ a/3 u Ac 7 ly g Q e /3 u^ , (8) 

M(E: -> S c7 ) = ier t3 e ta/aP u Tie Yq a ^i4 t , (9) 

where and g M are the photon polarization vector and momentum respectively, e is the 
charge unit. 

In order to derive the sum rules for the coupling constants we consider the correlator 
d A x e- ik - x {y(q)\T( V x(0)f} A (x)) |0> = e^ 7 f7 5 W^ a A ff G E , A ( W , W ') , (10) 



d'x e-^( 1 (q)\T( V ^(0)r ]A (x))\0) = 
^(e a q u - e u q a )^^(-gr + ^lht)^uaf3ae a q /3 v a G^ A (u,uj') , (11) 

d A x e-**< 7 ( g )|r(i^.(0)^( i r))|0> = 
e^— ltlh{-9t V + -^ltlt){e a q v - e v q a )Gz* , u') , (12) 

where k! = k — q, q 1 ^ = g M — (q ■ v)v fl , to = 2v ■ k, to' = 2v ■ k' and q 2 = 0. 

Let us first consider the functions G^ a (uj, oj') etc in (IIDp-(O). As functions of two 
variables, they have the following pole terms from double dispersion relation 

~ Al71lhh + ttA—, + 7rr~~ — , (13) 



(2A S - oj') (2A A - u) 2A S - 00' 2A A - uo 

Al7]2 + 7TT 7 + 7TT~— J (14) 



y/3 (2A2* — u')(2A A — lu) 2A s *-u' 2A A -u 
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4?/3 _ h*h_ c d 

y/3 (2A E * - u') (2A E - u) 2As* - w' 2A S - ' 1 ' 

where /s* etc are constants defined in As* = m^. — tuq. 

Keeping the two particle component of the photon wave function, the expression for 
G2* s a((jJ,u') with the tensor structure reads 



oc 



- 2 / dt J dxe~ lkx 5(-x - ^)Tr{[C(7(g)|M(0)M(x)|0) T C 7 ^S(-x)7 5 ] 

+ [CzS T (-x)C lll ]( 1 (q)\u(0)u(x)\0} l5 ]} , (16) 

where iS(—x) is the full light quark propagator with both perturbative term and contri- 
bution from vacuum fields. 



iS(x) = (0\T[q(x), g(0)]|0) = - M - . G g> 



"^-7^2 /' ^4°" • G M - 4m%G^( M x) 7 J + • • • • (17) 

l07T z JO 2 X A 

The light cone two-particle photon wave functions are H|: 



< 7(g)k(z)<V9(0)|0 >= ^e(gg) / due mqx {(e^q v - e v q lt )\x<p(u) + x A gi(u)] 

Jo 

+ [(qx)(e fl x l/ -e„x fl ) + (ex)(x fl q u -x u q fl )-x 2 (e fl q„-e l/ q^]g2(u)} , (18) 

< l(q)\q(xh^ 5 q(0)\0 >= f - e e e e v q p x° C du^ uqx ^{u) . (19) 

4 Jo 

Due to the choice of the gauge x^A^x) = 0, the path-ordered gauge factor 
Pexp (ig s Jq 1 dux^A^ux)) has been omitted. The <j){u),^}{v) is associated with the leading 
twist two photon wave function, while gi(u) and g2{u) are twist-4 PWFs. All these PWFs 
are normalized to unity, Jq 1 du f{u) = 1. 

Expressing (16) with the photon wave functions, we arrive at: 

/■oo rl . i , I 

G^ a {uj,u') = -{e u - e d ) / dt / due^-e lu -{(qq)[—x<P(u) 

JO jo 11 t 6 



+4r(0i(«) - 92{u))} + L^ u )t{{qq) + t-{qg s a ■ Gq))} + ■■■ . (20) 
7r z r 24 lb 



Similarly we have, 

G S)A (a;,u/) = G^ A {u,u') , (21) 
G s ., e (u,l/) = (e u + e d ) J dt J q due^^e^ {^^(«) 
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r {{m) + ^{qgsV-Gq))[x<t>{u)+t\ gi {u)-g 2 {u)}} + --- , (22) 
b lb 

where (qq) = — (225MeV) 3 , (qg s <7 ■ Gq) = ml(qq), = 0.8GeV 2 . For large Euclidean 
values of u and u' this integral is dominated by the region of small t, therefore it can be 
approximated by the first a few terms. 
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After Wick rotations and making double Borel transformation with the variables lo 
and uj 1 the single-pole terms in (0)-([y|) are eliminated. Subtracting the continuum con- 
tribution which is modeled by the dispersion integral in the region uj, uj' > uu c , we arrive 
at: 

VifefA = "6^( e « ~ e d )ae^^{x0(Mo)T 4 /3(f ) 

-4biW-g 2 N]T 2 /i(f) + ^W(i-^)} , (23) 

mh*fA = S(e u - e d )ae^ k {x0MT 4 / 3 (f ) 

-A\g 1 (u )-g 2 (u Q )]T 2 f i m+^mu )(l-<)} , (24) 



Vsh*h = iJ(e„ + e rf )e^^{^( Mo )T 3 / 2 (f ) 

-i(i-fl)[x0K)T/o(f)-|(^iM-^M)]} , (25) 

n k 

where f n {x) = 1 — e~ x XI tt is the factor used to subtract the continuum, uj c is the 

fc=0 

continuum threshold. u = r + T > ^ = j^ffi > ^2 ar e the Borel parameters a = 
— (27i) 2 (qq). We have used the Borel transformation formula: iS^e"^ = 5 (a — ^). 

Due to the heavy quark symmetry, As = As* and /e = /s* in the limit mg — > oo. So 
from fl23|) and (0) we have r? 2 = V^i- For the decays £*° -> S° 7 and £*++ -> we 
need make replacement (e u + e<f) — > 2ed, 2e u in (p5|). 

3 Radiative decays of heavy mesons 

We shall confine ourselves to the lowest lying three doublets and consider all possible 
radiative decay processes among them in the leading order of 1/ rriQ expansion. Denote the 
doublet (1+ 2+) with j e = 3/2 by (B 1: B* 2 ), the doublet (0+ 1+) with j e = 1/2 by (B' Q , B[) 
and the doublet (0~, 1") by {B, B*). 

The interpolating currents are given in [|l^] as 



A\ a = \ ~ A fa? - \l?V t ) q , (26) 



: 2 



4 v ' \ 1 3 
T T (-i) ( 2 



4+f = V 2 ~ K 2 ( 7rP " + 7?2Z>r ~ 3^ ia2Pt ) 9 ' (27) 
= /I ~ h ^tQ , 4^1 = I/I M*9 > (28) 



< + ,f - 771 M > j r + ,i = 775 ^Tf ? • (29) 
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(1+2+) -(0-1-)+ 7 
The decay amplitudes are 

M(B 1 - ZT 7 ) = eqee^e^aiie^q? + (a - /3M<?L(£?i, F") 

+[W*(« - - (« ~ + e mPa g s (B 1 , B*)} , (30) 

where the tensor structure associated with g l D (Bi,B*) and g 2 D (Bi,B*) is symmetric 
and antisymmetric under the exchange of (a. — 0) respectively. 

M(B 1 - B T ) = e ? ee^ a {(g«gf - \<$gf)g D {B x , B) + <?f<7s(£i, £)} , (31) 

M(B* - B 7 ) = e 9 ee^ aia2 g^[e^ Ql Y + (ai - a 2 )\g D {B* 2: B) , (32) 

M(B* - B* 7 ) = e^^ef^gf - ^g^ 2 )^^*, B*) + [e^(gfgf» 
-^tV as ) + («i - « 2 )]^(5 2 *, 5*) + [« 2/3 + (ai - a 2 )]<7s(£ 2 *, £*)} , 

where 77^, 7/ M and e M are polarization tensors for states 2 + , 1 + and 1" respectively 
and e q e is the light quark electric charge. 

Due to heavy quark symmetry, there exist only two independent coupling con- 
stants for the D-wave and S-wave decay respectively. Let ga = go{Bl,B) and 
g s = —g s {Bl,B*). Then we have: 

^g s {B u B*) = ^-g s {B u B) = g s , (33) 

^-gUB^B*) = V6g 2 D (Bi,B*) = ^-g D (B u B) 

= ~g 1 D (B*,B*)=gl(B*,B*)=g d . (34) 

The above relation is confirmed by our detailed calculation. 

For deriving the sum rules for the coupling constants we consider the correlator 

d\ e-^^T^ir^iCo)^^*)) |0> 

= e q e{e}(q?q? - \q 2 t gf)G D BlB {uj, J) + e* t a G s BlB (u;, J)} . (35) 

D S 

The functions G b ' iB (uj,uj ) in (gg) have the following double dispersion relation 
f- A U,%9Df{B x B) 

+ TTs ~, + ' > ( 36 ) 



(2A_ 1 -o/)(2A, 1 -w) 2A_i-cj' 2A , 3 - 

'2 """'a '5 
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where Apj e = mpj e — rriQ and fpj e are constants defined as: 

(0\J$j?{0)\j , ,P'dt) = fpA'SppStfTf*"-** . (37) 
Applying the same procedure as in section |2] we obtain 

+{qq)[ X <i>(u) + t 2 ( 9l (u) - g 2 (u))]} + ■■■ , (38) 
G s BiB (u,u') = - 2 -q 2 t G D BlB (u,u>) = l -{u-u'fG D BlB (u^') . (39) 

Finally we have: 

I A -,i +A +,3 a uj c A 

9df-,if+,\ = ge T {fuo^{u ) + —u [x(p(uo)Tf ( 1 ^)--(g 1 (u ) - g 2 (u ))]} , 



(40) 

d? ((ugijll) - Ug 2 (u))) ^ C , R| 

_4 J MjrJJI Im=«q 

Here we have used integration by parts to absorb the factor (q ■ v) 2 , which leads to 
the second derivatives in (f|l|). In this way we arrive at the simple form after double 
Borel transformation. 

(0+1+) -(0-,l-) + 7 

There exists only one independent coupling constant, corresponding to S-wave decay. 
The decay amplitudes are: 

M(B[ - 5* 7 ) = e q ee^e^ aV ' a e* g s (B[, B*) , (42) 

where rj' a is the polarization vector of B[. 

M(B[ - fl 7 ) = e q ee a V ' a g s (B[, B) , (43) 

M(B' - B*j) = e q eePepg s (B' , B*) . (44) 

The process B' — > B^ is forbidden due to parity and angular momentum conserva- 
tion. Due to heavy quark symmetry, we have 

g s (B[, B*) = g s (B[, B) = -g s (B' , B*) = 9l . (45) 

We consider the correlator 

d 4 x e -**< 7 (g)|T (j Q>+ i (0) J^i^)) |0) = e q ee?G %B * (cu, u') , (46) 
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where 

G B ' q b*{u,u') — --{qq)(q-v) jf dt ^ due^- u ^e^{ X( f>(u)+t 2 g 1 (u)} + . ■ ■ . (47) 
it a A --? +A +>? , d<p(u) ^ .. ,uj c , dgAu),. 

= r {x-sr^tT'- 4 * >!•-•■ (48) 

. (l + ,2+)^(0+,l+) + 7 

There exists only one independent coupling constant, corresponding to P-wave decay. 
The decay amplitudes are: 

+(<£gS' l -d9? l )9p(B 1 ,B' 1 )} , (49) 
M(Sx - 5 o7 ) = e,e Wa ?f e*^gV(£i, B ) , (50) 
M(5* - B o7 ) = e q ee^ a2 [{q^gf - jtfg?**) + (a, <-> a 2 )]g P (B* : B' ) , (51) 

M(B* 2 -> # l7 ) = e g ee^?/V^ QlQ J(g^r - ^« lQ2 ) + («i ~ « 2 )]<?p(5 2 *, *i) • 

(52) 

Due to heavy quark symmetry we have 
. fa 

^-gUBi, B[) = VQgUB^ B[) = v^p^i, B' ) = g P (B* 2 , B' ) = g P (B* 2 , B[) = g 2 . 

(53) 

We consider the correlator 

/ d'x e- ik -*(n(q)\T (0) J^7(x)) |0> = 

e q ee»l(qrgr - ^9^) + («i - « 2 )]G S|B , ( w , a/) , (54) 

where 
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1 /*00 /" 1 / 

= —<??>(?•«)/ due^V^) + tV («)} + ... . 

O JO JO 

(55) 

~32tt 2 iX du h{ T> du ilu=u °- [ ' 

M(B[ -> 5 o7 ) = e^ieiw . (57) 



In order to derive #3, we consider the correlator 

J d 4 x e-^( 7 (g)|T (j 0>+ ,i(0)^ +i i(x)) |0) = e q ee a ^e,q u v a G 3^(00^') , (58) 
where 

G BiB ,(u,u') = ll°°dt j'due^fe^ijMu) + (qq)[ X <f>(u) 

+t 2 ( 9l (u) - g 2 (u))}} + ■ ■ ■ . (59) 

1 2A 1 A 

9sfl h = --e^{-Mu ) + ^[X0K,)T/ O (^) - - ( 9l (u ) - g 2 (u ))]} . (60) 
• B* -> £7 

M(£* - £ 7 ) = e g ee a ^e a e;g^^ 4 . (61) 
In order to derive #4, we consider the correlator 

J d A x e- ifc -*( 7 (g)|T (^,-,1(0)^1^)) |0) = e q e^e^q u v a G B * B {u;,u') , (62) 

where 

GW^V) = 7 / A / ^"^^{--M^) + (gg)M(u) 
4 jo jo 4 

+i 2 teW- £2 («))]} + ■•• • (63) 

Note this coupling was calculated in H using LCQSR. But there the contribution 
from the photon wave function ip(u) has not been taken into account. 

1 2A _ 1 A 

94fl 1 = —e^WM + ^[X0K)T/ O (^) - = G&(«o) - fls(«o))]} ■ (64) 

'2 4 Z7T Z i i 

. B* - B l7 

M(B* - B l7 ) = e ? ee^% 9 ^ W2 ^ • 

(65) 

In order to derive g 5 , we consider the correlator 

/ d*x e -*»< 7 (?)|T (^ +f} (0)J^|(ar)) |0) = 
e g ee^%^(sr<zr + <?r? t ai - ^rV)(2gf<T + gf \ P t)G b*b^ > «0 , (66) 
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where 

G B * Bl (u,u') = r dt f'due^-^e^uil - u ){~fi/>(u) 

16 Jo Jo 4 

+(qq) \x<P(u) + t 2 ( 9l (u) - g 2 (u) )]} + •• • . (67) 



/g 2A + 3 <x ui \ 

9sf + ,l = -j^e^-uoil - u ){f4j(u ) + —[ X (f){uo)Tf (^) - — (gx(u ) -02M)]} . 

(68) 

4 Determination of the parameters 

The leading photon wave functions receive only small corrections from the higher conformal 
spins || so they do not deviate much from the asymptotic form. We shall use [|J 

4>{u) = 6uu , (69) 

= 1 , (70) 

9i(u) = ~u(3-u) , (71) 

92{u) = -\u 2 . (72) 

with / = 0.028GeV 2 and \ = -4.4GeV 2 pf at the scale fi = lGeV. Using this value of x, 
the octet, decuplet and heavy baryon magnetic moments have been calculated to a good 
accuracy p4] p|, p^]. 

We need the mass parameters A's and the coupling constants /'s of the corresponding 



interpolating currents in the leading order of a s as input. The results are [JTT 

A A = 0.8 GeV f A = (0.018 ± 0.002) GeV 3 , 

A s = 1.0 GeV / s = (0.04 ± 0.004) GeV 3 . (73) 

A +)3/2 = 0.82 GeV / +i3/2 = 0.19 ± 0.03 GeV 5/2 , 
A+,1/2 = 1.1 GeV / +! i /2 = 0.40 ± 0.06 GeV 3/2 , 

A_ 1/2 = 0.5 GeV /_ i /2 = 0.25 GeV 3/2 . (74) 

We choose to work at the symmetric point T\ = T 2 = 2T, i.e., uq — |. Such a choice 
is very reasonable for the symmetric sum rules (|25|) , (|60|), (|64D and ( |68D since S* and E c , 
and the three meson doublets are degenerate in the leading order of HQET. Moreover, the 
mass difference between E* and A c is only about 0.2GeV. The (0 + , 1 + ) doublet lies only 
slightly below (1 + , 2 + ) doublet. Due to the large values of T\, T 2 used below, the choice of 
T\ = T 2 is also reasonable for sum rules ( p3| ) and (|56|). 

Note the choice T\ = T 2 is not unique for the asymmetric sum rules (0), (|4l|) and (|48D 
since the initial and final mesons have different masses. But the choice T\ = T 2 will enable 
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the clean subtraction of the continuum contribution, which is cruicial for the numerical 
analysis of the sum rules. In our case the sum rules are stable with reasonable variations 
of the Borel parameter 7\ and T 2 . Such a choice does not alter significantly the numerical 
results. Based on these considerations we adopt w = \ for for the sum rules fl40|), (f4~l|) 
and (P|) too. 



5 Numerical results and discussion 

5.1 Numerical analysis of the baryon sum rules 

We now turn to the numerical evaluation of the sum rules for the coupling constants. 
Since the spectral density of the sum rule (p3|)-(p5|) p(s) is either proptional to s 2 or s 3 , the 
continuum has to be subtracted carefully. We use the value of the continuum threshold 
u c determined from the corresponding mass sum rule at the leading order of a s and 1/tbq 

III- 

The lower limit of T is determined by the requirement that the terms of higher twists 
in the operator expansion is reasonably smaller than the leading twist, say < 1/3 of the 
latter. This leads to T > 1.3 GeV for the sum rules (P3"|)-([Z5]). In fact the twist-four terms 
contribute only a few percent to the sum rules. The upper limit of T is constrained by 
the requirement that the continuum contribution is less than 50%. This corresponds to 
T < 2.2GeV. 

The variation of 771,3 with the Borel parameter T and u c is presented in FIG. 1 and 
FIG. 2. The curves correspond to uo c = 2.4, 2.5, 2.6GeV from bottom to top respectively. 
Stability develops for the sum rules (|23"D and (|25|) in the region 1.3 GeV <T<2.2 GeV, we 
get: 

VifefA = (7.0 ± 0.9) x 10~ 4 GeV 5 , (75) 
Vsfafc = (3.9 ± 0.5) x 10- 4 GeV 5 , (76) 

where the errors refers to the variations with T and uj c in this region. And the central 
value corresponds to T = 1.6GeV and u c = 2.5GeV. 
Combining ([73]) we arrive at 

= (1.0±0.2)GeV _1 , (77) 
773 = (0.24 ± 0.05)GeV _1 . (78) 

5.2 Numerical analysis of the meson sum rules 

We now turn to the numerical evaluation of the sum rules for the coupling constants. The 
lower limit of T is determined by the requirement that the terms of higher twists in the 
operator expansion is less than one third of the whole sum rule. This leads to T > 1.0 
GeV for the sum rules (@g), f|), (0), (§4]) and In fact the twist-four 

terms contribute only a few percent to the sum rules for such T values. The upper limit 
of T is constrained by the requirement that the continuum contribution is less than 30%. 
This corresponds to T < 2.5GeV. With the values of photon wave functions at uo = | we 
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obtain the left hand side of these sum rules as functions of T. The continuum threshold 
is uj c = 3.0 ± 0.2GeV except u c = 2.4 ± 0.2GeV for the sum rule (|64"D . Stability develops 
for the sum rules in the region 1.0 GeV <T<2.5 GeV. The results are shown in FIG. 3-9. 
Numerically we have: 





s = -(3.0 ±0.2) x 10~ 2 GeV 2 , 


(79) 


9sf-,if+ 


a = -(1.9 ±0.2) x 10" 2 GeV 4 , 

2 v ' 


(80) 


9if-,±f+ 


i = -(1.5 ±0.5) x 10~ 2 GeV 3 , 


(81) 


92f+lf+ 


3 = -(5.5 ±0.4) x 10~ 2 GeV 3 , 


(82) 


9*fl,\ = 


(0.28 ± 0.04) GeV 2 , 


(83) 


9,fl I = 

' 2 


(8.9 ± 0.5) x 10~ 2 GeV 2 , 


(84) 


9bfl 2 = 

^'2 


-(2.3 ±0.3) x 10" 2 GeV 2 , 


(85) 



where the errors refer to the variations with T in this region and the uncertainty in uo c . 
And the central value corresponds to T = 1.5 GeV and u c = 3.0 ± 0.2 GeV except that we 
use oj c = 2.4 ± 0.2GeV for the sum rule (H). 

With the central values of f 's in ( [74] ) we get the absolute value of the coupling constants: 

9d = -(0.63 ±0.10) GeV~ 2 , (86) 

g s = -(0.40 ± 0.05) , (87) 

Si = -(0.20 ±0.06) , (88) 

g 2 = -(0.72 ± 0.07) GeV" 1 , (89) 

g 3 = (1.8 ±0.3) GeV" 1 , (90) 

04 = (1-4 ± 0.2) GeV- 1 , (91) 

g 5 = -(0.64 ± 0.08) GeV" 3 . (92) 

Note we have only considered the uncertainty due to the variations of the Borel pa- 
rameter and the continuum threshold in the above expressions. There are other sources 
of uncertainty. The input parameters \ an d / are associated with the photon distribution 



amplitude. Especially the value of x has been estimated with QCD sum rules [|1^] and 
with the octet baryon magnetic moments as inputs using the external field method fll5 |. 
Both approaches yield consistent results x ~ —4.4 GeV. With this value the octet, decu- 
plet and heavy baryon magnetic moments derived using the external field method are in 
good agreement with the experimental data. So we expect its accuracy is better than 30%. 
The value of / has been estimated with the vector meson dominance model also with an 
accuracy of 30% [0]. 

The light cone sum rules for the coupling constants gi and the mass sum rules for the 
heavy hadrons in HQET both receive large perturbative QCD corrections. But their ratio 



does not depend on radiative corrections strongly because of large cancellation [|17]] . In the 
present case, the uncertainty of the coupling constants gi due to radiative corrections is 
expected to around 10% while the couplings f\ etc are affected significantly. 

Another possible source of error is the truncation of the light cone expansion at twist 
four. We take the sum rules for r/i for example. At T = 1.5 GeV, the twist-four term 
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involved with gi,g% is only —2.5% of the leading twist term after making double Borel 
transformation to (p0|). Even after the subtraction of the continuum and excited states 
contribution the twist-four term is only —15% of the twist-two one in (|2"5|). So the light 
cone expansion converges quickly. We expect the contribution of higher twist distribution 
amplitudes to be small. 

We have calculated the coupling constant in the leading order of HQET. The 1/toq 
correction is sizable for the charm system. But for the bottom system the 1 / mq correction 
is typically around 5% ~ 10% for the pionic coupling constants ||. We expect the 1/mq 
correction to the electromagnetic coupling constants is of the same order. The inherent 
uncertainty of the method of QCD sum rules is not included, which is typically about 10%. 

5.3 Decay widths of heavy hadrons 

With these coupling constants we can calculate the decay widths of heavy hadrons. 
The decay width formulas in the leading order of HQET are 



r(s fe - 


-> A h7 ) = 


4?7i«|gl 3 , 


r(s*- 


-> Afc7) = 


2 3 3m 2 + m} 
^ 3mf ' 


nn- 


- s 67 ) = 


3m?+roJ 

mm 3m 2 » 


r(£?i- 


- B*l) = 


2 14 


r(£?i- 


-^7) = 


Y e l a (j^9 2 d \q\ 5 + g 2 M) , 


t(b; - 


-^7) = 


2 

^ q ag d \q\ , 


r(B* - 


- B*i) = 


~- e 2 q a(^g 2 d \q\ 5 + Ag 2 s \q\) , 


T(B[- 


- B*l) = 




T(B[- 


-57) = 






- B*l) = 


3 9 2 1 — h 


T(B* - 


-57) = 


2 2 1 — h S 



where \q\ = l 2m 1 , m^ mj is the parent and decay heavy hadron mass. 

We apply the leading order formulas obtained above to the excited states of bottomed 
hadrons using the central values of the coupling constants in the previous section. 

r(E^A t7 ) = 1 31x( I J|-^) 3 keV, 
r( E^A t7 ) = 313xMM\eV, 
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r(Er Kl) = 2-2 x (^i^) 3 keV , 
r(E^EM=0.14x( s iL)\eV, 
r(EJ" -> S^ 7 ) = 0.56 x ( — i|L - j keV 



63.4MeV 



r ( zr^^) = i.4x^^kev, 

IY£° -> fi° 7 ) = 84.4 x f — S — ] keV , 

F(B+ -> 5+ 7 ) = 338 x f l|L— : ] keV , 

v 1 /; V490MeV/ 

-> 5*° 7 ) = 42.2 x [ — JfL— ] keV , 
v 1 /; V377MeVy 

r(5+ -> 5*+ 7 ) = 169 x [ — M— - ] keV , 
v 1 /; V377MeV; 

T(B;° -> fi° 7 ) =5.8x1 S ] keV , 

V 2 [; v 537MeV / 

r(5 2 * + -> = 23 x ( ^ ] keV , 

V 2 [; v 537MeV / 

r(B- - B '+ 7 ) = 844 x (jjL) keV . (94) 

The uncertainty of the decay width is typically about 30%. 

We do not present numerical results for the radiative decay widths for the charmed 
hadrons since 1/mq corrections are sizable for the charm system while such corrections are 
only a few percent of the leading order term for the bottom system ||. 

In summary we have calculated the coupling constants of photons with the heavy 
baryons and the lowest three heavy meson doublets using the light cone QCD sum rules 
with the photon wave functions in the leading order of HQET. We hope these calculations 
will be tested in the future experiments. 
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Figure Captions 



FIG. 1. Dependence of fsfwi on t ne Borel parameter T for different 
values of the continuum threshold lo c . From top to bottom the curves 
correspond to u> c = 2.6, 2.5, 2.4 GeV. 
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